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We present polynomial Poisson algebras for the 8 classical potentials in 
two-dimensional Euclidian space that separate in cartesian coordinates and 
allow a third order integral of motion. Some of the classical superintegrable 
potentials do not coincide with quantum ones, but are their singular limits. 
We present the trajectories for all these classical potentials. We find that all 
bounded trajectories are periodic. 
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1 Introduction 



The purpose of this article is to study the properties of a certain class of 
superintegrable systems in a real two-dimensional Euclidian space E 2 . These 
are classical hamiltonian systems that allow three functionally independent 
integrals of motion. One of them is the Hamiltonian H. The other two are 
second and third order polynomials in the momenta, respectively. Moreover, 
in most of the article we shall restrict to the case when the potential allows 
separation of variables in Cartesian coordinates, so that we have 

H=\{Pl + P 2 2 ) + f{x)+g{y) (1.1) 
A=\(P?-P2) + f(x)-9(v) (1-2) 

and X 3 will have the general form 

^3= ^jkLlP{P 2 k + k 1 (x,y)P l + k 2 (x,y)P 2 , (1.3) 

i+j+k=3 

where are constant and L 3 is the angular momentum 

L 3 = xP 2 - yP x . (1.4) 

Superintegrable systems are Hamiltonian systems with more integrals of 
motion than degrees of freedom. In dimension n=2 the maximal possible 
number of functionally independent integrals is 2n-l=3 and that is the case 
we are considering. 

A large body of literature exists on superintegrable systems, inspired by the 
superintegrability of the harmonic oscillator 1 ' 2 and Kepler-Coulomb 
potentials 3 ' 4 . A systematic search for superintegrable systems in E 2 and E 3 
was conducted in the 1960ties 5,6 ' 7 . Like much of the later work on 
superintegrable systems it was restricted to the case of second order 
integrals of motion 8 ' -' 13 . This case turned out to have an intimate 
connection with the separation of variables in the Hamilton- Jacobi and 
Schrodinger equations. 

An interesting discovery was that in general the quadratic integrals of 
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motion do not generate finite dimensional Lie algebras, but more 
complicated algebraic structures, namely quadratic algebras 14-18 . 
Much less is known about integrable and superintegrable systems with 
third and higher order integrals of motion. In 1935 Drach found 10 different 
potentials in a complex Euclidian plane allowing a third order integral of 
motion in classical mechanics 19 ' 20 . In 1984 Hietarinta 21 showed that in the 
case of integrals that are third, or higher order polynomials in momenta, 
quantum and classically integrable potentials may not coincide. 
A systematic search for superintegrable systems involving at least one third 
order integral was started in 2002 22 . It was shown that in the case of one 
first order and one third order integral the only classical superintegrable 
systems were known ones (like V = ar 2 or V — or -1 ) for which the third 
order integral was the product of a first and second order one. However, in 
the quantum case a new superintegrable potential of this type exists, 
namely Vi(x,y) = (hu) 2 k 2 sn 2 (ujx, k) where sn(ux, k) is a Jacobi elliptic 
function. The classical limit (H — > 0) is free motion ! 22 
All classical and all quantum potentials allowing a second order integral of 
the form (1.2) and a third order integral (of any form ) were found in ref. 23 . 
The classical ones are all expressed in terms of elementary functions. The 
quantum ones are quite different and some of them involve elliptic functions 
and Painleve transcendents 24 . We leave the quantum case for a future study 
and concentrate here on the classical superintegrable Hamiltonian systems 
of the form (1.1),.. .,(1.3) 

In Section 2 we consider a general superintegrable Hamiltonian system in a 
two-dimensional space (not necessarily Euclidian) with a scalar potential. It 
allows one quadratic and one cubic integral of motion. We construct the 
most general cubic Poisson algebra generated by these integrals. We express 
the Casimir operator K of this algebra in terms of the Hamiltonian H. The 
general formalism of Section 2 is applied in Section 3 to the 8 existing 
superintegrable systems in E 2 allowing separation of variables in Cartesian 
coordinates and having an additional third order integral of motion. 
Solutions of the equation of motion and trajectories in these 8 systems are 
discussed in Section 4. Section 5 is devoted to conclusions. 
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2 Cubic Poisson algebras 

We start out with a more general Hamiltonian system than (1.1)... (1.3). We 
will not assume that we have a natural Hamiltonian in E 2 , but that we 
have a superintegrable system with a quadratic Hamiltonian and one 
second order and one third order integral of motion. We put 



H = a(q u q 2 )P 2 + 2b(q u g 2 )PiP 2 + c(q u q 2 )P 2 + V(q u q 2 ) 

A = A( qi , q 2 , P u P 2 ) = d(q u q 2 )P 2 + 2e(q u q 2 )PiP 2 + 
f(qi, l2)Pl + g(qi, q 2 )Pi + h( qi , q 2 )P 2 + Q(q u q 2 ) (2.1) 
B = B( gi , q 2 , P u P 2 ) = u(q u q 2 )P? + 3v(q u q 2 )P?P 2 + 3w( qi , g 2 )PiP| 
+x(q u q 2 )P^+j(qi, g 2 )P 1 2 +2/c(g 1 , q 2 )PiP 2 +l(qi, q 2 )P 2 +m( qi , q 2 )Pi+n(q 1 , q 2 )P 2 +s( qi , q 2 ) 



where p and qi are canonical momenta and coordinates and 

functions to be determined. Since A and B are integrals of motion we have. 



where {, } is the Poisson bracket 

Since B is cubic in the momenta, we cannot expect {H, A and B} to 
generate a quadratic algebra, like those obtained in 14 ' - ' 18 . We will however 
try to close the algebra at the lowest order possible, namely 3. We put 



{H,A} = {H,B} = 



(2.2) 



{A,B} = C 
{A, C} = aA 2 + 2(3 AB + ^A + 5B + e 
{B, C} = fiA 3 + vA 2 + pB 2 + 2a AB + £A + r]B + ( 



(2.3) 



Since we have {C, {A, P}}=0 the Jacobi identity reduces to 



{A,{B,C}} = {B,{A,C}} 



(2.4) 



The Jacobi identity implies 
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p = —(3, a = —a, r] = —7 and we obtain the cubic algebra. 

{A,B} = C 

{A,C} = aA 2 + 2(3AB + 1 A + 5B + e (2.5) 
{B, C} = fiA 3 + vA 2 - (3B 2 - 2aAB + £A - 1 B + (. 

The coefficients a , j3 and \i are constants, but the other ones can be 
polynomials in the Hamiltonian H. The degrees of these polynomials are 
dictated by the fact that H and A are second order polynomials in the 
momenta and B is a third order one. Hence C can be a fourth order 
polynomial. We have 

a = a , (3 = /3o, \i = Mo (2.6) 
7 = 7o + 7i#, 5 = 5 + 5iH, e = e + e x H + e 2 H 2 
u = u Q + u l H,i = i + i 1 H + i 2 H 2 
( = (o + ( 1 H + ( 2 H 2 + ( 3 H 3 , 

where ao,...,C3 are constants. A Casimir operator K of a polynomial algebra 
is defined as an operator Poisson commuting with all elements of the 
algebra. For the algebra (2.5) this means 

{K,A} = {K,B} = {K,C} = (2.7) 

and this implies 

K = C 2 -2aA 2 B-2f3AB 2 -2-fAB-5B 2 -2eB + -fiA 4 + -vA 3 + £A 2 + 2(A. 

(2.8) 

Thus K is a polynomial of order 8 in the momenta, since the Hamiltonian H 
also satisfies relations (2.7) we can expect K to be a polynomial in H and 
we write 

K = k + k 1 H + k 2 H 2 + k 3 H 3 + k 4 H* (2.9) 
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where ko, .., k± are constants. Note that eq. (2.9) together with (2.8) 
represents a polynomial relation between the integrals H, A, B and C in 
agreement with the fact that only 3 of them can be functionally 
independent. 

3 Cubic Poisson algebras for classical 
superintegrable systems 

Let us now apply the formalism of Section 2 to all classical superintegrable 

potentials that separate in Cartesian coordinates and allow a third order 

integral of motion 22 ' 23 . They have the form (1.1),.. .,(1.3) , so we have 

a (<?i> Q2) = c(qi, q 2 ) = 1, b(qi, 52) = in (2.1). Eight such systems exist, four 

of them are well-known, four first presented in 23 . 

Case 1. The isotropic harmonic oscillator. 

We have 



The formalism of Section 2 applies, but we cannot expect it to provide 
anything new. indeed, in this case we have a well known u(2) algebra of 
first and second order integrals of motion : 




(3.1) 



L 3 = xP 2 - yP u X 1 = P 2 + u 2 x 2 , X 2 = P 2 2 + u 2 y 2 
X 3 = P X P 2 + oo 2 xy 



(3.2) 



we can chose 



A = P 2 - P| + co 2 (x 2 - y 2 ) =X 1 -X 2 
B = L 3 X 2 = xP 3 - yP x P 2 - uJ 2 y*P 1 + u 2 xy 2 P 2 



(3.3) 



we obtain 



{A, B} = C = -4X 2 X 3 
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{A,C} = -16u 2 B 
{B, C} = 2A 3 - QHA 2 + 8H 3 



(3.4) 



Thus, relations (2.5) are satisfied, but we remain in the enveloping algebra 
of u(2). The Casimir operator of the cubic algebra in this case is 

K = C 2 + lQu 2 B 2 + 2A A - 6HA 3 + 8H 3 A = 16H 4 (3.5) 

However, in this case H is simply the central element of the u(2) algebra 
(3.2) and it generates the center of the enveloping algebra of u(2)( a hence 
also of the cubic algebra) . 



Case 2. A quadratically superintegrable Hamiltonian 



5,6 



The Hamiltonian allows two second order integrals of motion, namely 



X 1 = P 2 -P 2 + u 2 (x 2 -y 2 ) + ^- 2 ^ (3.7) 

x y 



X 2 = L 2 3 + 2r 2 (- + -). 



b c 

~x 2 y 2 

The Hamiltonian H allows separation of variables in cartesian and polar 
coordinates and also in elliptic ones. We shall identify X± — A and the 
cubic integral B is equal to 



— +u;x)P 2 -xy( — 



B = ~{X 1 ,X 2 } = xP 1 P 2 2 - yP 2 P 2 + xy(— + uo 2 x)P 2 - xy{— + oo 2 y)P l 

(3.8) 



{A, B} = C = A(A 2 - 2HA - H 2 + 8oo 2 X 2 - \Quo 2 (b + c) 

{A, C} = -QAlo 2 B (3.9) 
{B, C} = -2A 3 + 8H 2 A + 64(c - b)u 2 H + 32(c + b)u 2 A 
K = C 2 + 6Acu 2 B 2 - A 4 + (8H 2 + 32(c + b)u 2 )A 2 + 128(c - b)u 2 HA (3.10) 
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= 16H 4 - 128u 2 (c + b)H 2 + 1024c/&c 



In this case there is no underlying finite-dimensional Lie algebra, however 
Xi,X 2 and H generate a quadratic Lie algebra 16,18 

{Xi, X2} = C q 

{X 1} C q } = -8X 2 + 16Xi# - \Qu 2 X 2 + 32u 2 {b + c) (3.11) 
{X 2 , C q } = 16X 1 X 2 - 16HX 2 + 32H(c - b) 

With C q = —B and the relations (3.9) of the cubic algebra are 
consequences of (3.11) 

Case 3. A further quadratically superintegrable Hamiltonian 5 ' 6 is 

H=S- + ?l + ^(4x 2 + y 2 ) + ± + cx (3.12) 

The two second order integrals of motion are 

A = Pi - Pi + uj 2 (Ax 2 - y 2 ) - p + 2cx (3.13) 

X 2 = xPl - yP x P 2 + ^ - 2cm/ 2 - °-y 2 

y l 2 

and the Hamiltonian allows the separation of variables in cartesian and 
parabolic coordinates. The third order integral is 

B = X 2 } = P x Pl + (^ - u; 2 y 2 )P 1 + (Aoo 2 xy + cy)P 2 (3.14) 

4 y z 



{A, B} = C = l6u 2 X 2 
{A, C} = -6Alo 2 B (3.15) 
{B, C} = 1Qlo 2 HA + 1Qlo 2 H 2 - \2u 2 A 2 + 8c 2 H - 4c 2 A + 128c/6 
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K = C 2 + 64u 2 B 2 -8u 2 A 3 + ( y 16u 2 H-4c 2 )A 2 + ( y 32uj 2 H 2 + 16c 2 H + 256uj 4 b)A 

(3.16) 

= 64u 2 H 3 + 16c 2 H 2 - 512u 4 bH - 128u 2 bc 2 



As in Case 2 we have a cubic algebra but all relations in it follow from the 
already known quadratic algebra 16,18 . 

Case 4. We have 

H=S- + 5l + ^(9x 2 + y 2 ) (3.17) 

Which is anisotropic harmonic oscillator with a rational frequency ratio 
cji/c<j 2 =3/1. The second and third order integrals are 

A = P 2 -P 2 + u 2 (9x 2 -y 2 ) (3.18) 
B = -yPxPl + xPl + \u 2 y z Pi - 3u 2 xy 2 P 2 



{A, B} = C = -4PiP| + V2u 2 y 2 P x P 2 - 36u 2 xyP% + 12u 4 xy 3 ) 

{A,C} = -UWB (3.19) 
{B, C} = 2A 3 + 8H 3 - QHA 2 

K = C 2 + UAlu 2 B 2 + A 4 - AHA 3 + 1QAH 3 = 1QH 4 (3.20) 

In this case the integral B is irreducible, i.e. it is not the Poisson bracket of 
lower order integrals. The cubic algebra (3.20) is related to a u(2) 
invariance algebra constructed by Jauch and Hill 1 for any anisotropic 
harmonic oscillator in E 2 with a rational frequency ratio. Their Lie algebra 
for the ratio Ui/u2=3/l is generated by the integrals. 

{H, B i = J^A' B2 =B~A> Bs = A} (3 - 21) 

It has actually been shown that in many cases in classical mechanics 
functions of integrals of motion can be constructed that generate finite 



9 



dimensional Lie algebras 25 . In the case of the anisotropic oscillator (3.17) it 
suffices to take the fractions B\ and B 2 of the polynomials (3.17),.. .,(3. 20). 
The Hamiltonians considered so far namely (3.1), (3. 6), (3. 12) and (3.17) are 
all superintegrable, both in the classical and quantum cases 23 . 

The remaining 4 cases are quite different in that the systems are integrable 
only in the classical case. They are all obtained as singular limits of 
quantum integrable systems. As we shall see, the quantum systems are 
quite different and the potentials are expressed in terms of Painleve 
transcendents. 



Case 5. The classical Hamiltonian and two integrals of motion in this case 
are 



h = S. + ^+PlVW\+PiVW\,A = ^--^ +p 2 1 VW\-PiVW\ (3-22) 



B = /3 2 4 P 3 + e/3 4 P| + 3/%PW\x\Pi + eSPtPW\v\P2 



?4 D 3 



e=l,xy>0 (3.23) 
e = — 1, xy < . 

In the quantum case an integral of the type B exists if the potential 

V(x,y)=\/i(:r) + V 2 (y) satisfies 23 

h 2 Vl'{x) = W?(x) - 6ftx (3.24) 

The classical (and singular) limit h — > yields the potential in (3.22). The 
cubic algebra (2.5) in this case simplifies and we get 

{A, B} = C = 3(3tti 

{A,C} = (3.25) 
{B,C} = 
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K = 9(3^1 



(3.26) 



Thus {A,B,C} actually generate a nilpotent Lie algebra, isomorphic to the 
Heisenberg algebra in one dimension (since C is constant). 

We can see from eq. (3.24) that the quantum case will be completely 
different. Indeed eq. (3.24) can be solved in terms of the first Painleve 
transcendent P 2i . 

Case 6. The classical Hamiltonian and two integrals of motion are 

p2 p2 , ,2 p2 p2 , ,2 

H = f + f + yy 2 + V(x), A = ff-ff-^ + V(x) (3.27) 



B = -yPf + xP 2 P 2 + (^-x 2 - W)yP l - ~ 2 (^x 2 - 3V)V X P 2 
where V satisfies a quartic equation 

-9V(x) 4 +Uuj 2 x 2 V(x) 3 + (Qd-3^-x 4 )V(x) 2 + (^x 6 -2io 2 x 2 )V(x) (3.28) 

+ (cx 2 -d-d l 4x*-^:x 8 ) = 
2 lb 

In the quantum case V satifies a fourth order differential equation 23 

frV (4) = 12u 2 xV + 6(V 2 )" - 2uj 2 x 2 V" + 2u 4 x 2 (3.29) 

that can be solved in terms of the Painleve transcendent Piv 2A - Eq. (3.27) 
is the solution of (3.29) for h — > and c and d are integration constants. In 
general, eq. (3.27) has 4 roots and the expressions for them are quite 
complicated. A special case occurs if u; 2 ,c and d satisfy. 

2 3 LU 8 b 3 , oo 4 b 2 , x 

c = ^^> d = ^r ( 3 - 3 °) 

where b is an arbitrary constant. Then eq. (3.28) has a double root and we 
obtain 

co 2 



V lj2 (x) = —(2b + hx 2 ± Axy/b + x 2 ) (3.31) 
18 
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V s (x) = V 4 (x) = (-^ + ^x 2 ) (3.32) 
For V(x) satisfying (3.28) the cubic algebra is 

{A,B} = C 
{A, C} = -Auj 2 B (3.33) 

AV2. ,4 A L3, ,6 

{B, C} = 8A 3 + 12HA 2 - AH 3 - A^-A + 



27 729 

if = C 2 + cu 2 5 2 + AA A + - 4k4 2 + {-8H 3 + ^ )A (3.34) 

27 729 



Case 7. The classical Hamiltonian and two integrals of motion in this case 
are 

H = 3 + 3. + a 2 \y\ + b 2 ^\x],A = 3- - 3- - a 2 \y\ + b 2 ^/R (3.35) 
B = P 3 + 36 2 v^Pi + 

2cr 



e = -l,xy>0 (3.36) 
e = 1, xy < . 

In the quantum case the potential is V(x,y)=ay + V(x) where V(x) satisfies 

h 2 V" = 6V 2 - 6b 4 x, b ^ (3.37) 

This is solved in terms of the first Painleve transcendent for K ^ 0. For 
h — the singular limit is Y=b^. The cubic algebra in this case is very 
simple and reduces to a Lie algebra, the Heisenberg algebra as in Case 5. 
We have 

{A, B} = C = 3b 4 
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{A,C} = (3.38) 
{B,C} = 
K = 9b 8 (3.39) 



Case 8. In this case we have 



H=^- + ^- + ay + V(x), A = 3- - 3 - ay + V(x) (3.40) 
B = aPf - bP 2 P 2 + a(3V(x) - bx)P l - 2bV(x)P 2 



where V(x) satisfies a cubic equation 

V(x) 3 - 2bxV(x) 2 + b 2 x 4 V(x) - d = (3.41) 

This is again obtained from a singular limit of a second order nonlinear 
ODE in the quantum case. It is, for K ^ 0, solved in terms of the Painleve 
transcendent Pjj. The cubic algebra in this case reduces to 

{A,B} = C = 2ab(A + H) 

{A,C} = (3.42) 
{B,C} = —Aa 2 b 2 (A + H) 

K = Aa 2 b 2 H 2 (3.43) 

We see that (3.42) is actually a Lie algebra with centre H. If we put 

H+A 1 ^ , 

A l = —,B^--B (3.44) 

We obtain the solvable decomposable Lie algebra 

{H u B,} = B u {H u H} = {B 1 ,H} = (3.45) 
The cubic equation (3.41) can be taken to standard form by putting 

V(x) = y(x) + ^ (3.46) 
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We obtain 

y 3 + 3py + 2q = 0, p = -(y ) 2 , q = (y) 3 - ^ (3.47) 
This discriminant is 

. o q b doc , , 

^ = g 2 +p 3 = T -^ (3.48) 

For D < 0,i.e. x 3 > |^ with d<0 or x 3 < |^ with d>0 it has 3 real roots, 
For D > 0,i.e. x 3 < ||f with d<0 or x 3 > |pf with d>0 it has 1 real root yi 
and 2 complex ones,y 3 = £/ 2 . 

The corresponding potentials are 

, . 2te 2 1 / 3 &V 
Vi(a;) = + = - 3.49 



(27d - 2b 3 x 3 + 3 v / 3v / 27rf 2 - 46W) 1 / 3 
+ 32 1 / 3 
Tr/ s 2bx (l + iV3)b 2 x 2 

V 2 (x) = + ^ ' - 3.50 

3 32 2 /3(2 + d- 26 3 x 3 + 3 v / 3V27d 2 - 46 3 dx 3 ) 1 /3 v ' 

(1 - iV?>){21d - 2fr 3 x 3 + 3y^3V27rf 2 - 46W) 1 / 3 
+ 62V 3 

, N 2bx (1 - i-v/^&V 

V 3 (x) = + ; ' ; 3.51 

V ; 3 32 2 / 3 (2 + d - 26 3 x 3 + 3VsV^7d 2 - 46 3 rfx 3 ) 1 /3 v ' 

(1 + iV3)(27d - 2& 3 a: 3 + 3y / 3y / 27d 2 - 4&W) 1 / 3 
+ 62V 3 

Multiple roots occur for D=0.This happens (for x ^ const) for d=0 only 
and then we have 

V = bx (3.52) 

as a double root (and V=0 as the simple one). We can consider the 
particular Hamiltonian, 

H=^ + ^ + a 2 \y\ + b 2 \x\,A = ^-^-a 2 \y\ + b 2 \x\ (3.53) 
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B = a 2 Pf - b 2 P 2 P 2 + 2a 2 6 2 |:r|P 1 - 2b 4 \x\P 2 



4 Trajectories for classical superintegrable 
systems 

If a Hamiltonian system is maximally superintegrable and satisfies certain 
analyticity properties, then all of its bounded trajectories are closed and 
the motion is periodic 26 . Here we shall discuss the trajectories for all the 8 
systems of Section 3. The trajectories can be obtained in a uniform manner 
directly from the integrals of motion. Indeed in all cases we can put 

\p?+m = E 1 (4.i) 

^Pi + g(y) = E 2 (4.2) 

B = /iPf + uP 2 P 2 + pP x Pl + aPl + 0P + tpP 2 = k . (4.3) 

In (4.3) jj, ,v ,p and a are low order polynomials in x and y, and ip are 
functions of x and y (all of them known). The constants E\ and E 2 are 
positive, k arbitrary. From (4.1) and (4.2) we obtain P 2 and P 2 in term of x 
and y, respectively. From eq(4.3) we obtain 

(P 2 (/iP 2 + pPl + 0) 2 - P 2 >P 2 + oPl + ij) 2 ) 2 + k 4 (4.4) 

= 2k 2 (P 2 (pP 2 + pP 2 + 0) 2 + P 2 2 {yPl + aPl + i,) 2 ) . 

Substituting for P 2 and P 2 from (4.1) and (4.2) we obtain the equation for 
the trajectories. Directly from eq (4.1) and (4.2) we see that the motion is 
bounded if there exist two constants, x Q and |/ ,such that we have 

f(x)>0, x 2 >xl g(y)>0,y 2 >y o . (4.5) 

This direct method of computing trajectories is universal, hence not 
necessarily convenient in any special case. For instance for the harmonic 
oscillator it gives a polynomial of order 16 in x and y which can then be 
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simplified to a second order one. However it avoids the problem of 
integrating the equation of motion and then eliminating time. 



Case 1. The trajectories of the harmonic oscillator are well known and are 
ellipses. 

Case 2. The trajectories are given in 6 and are easily obtained by 
integrating the equation of motion : 



E 1 E 2 2b . 



W = ~i + \l ~ —sin(2ut + Cl ) (4.6) 



CO 2 V 200^ UJ 



, „ E 2 El 2c . , , 

V(t) 2 = J + \I^I- - 2 sin(2ut + c 2 ) (4.7) 

Case 3. The trajectories are also given in 6 and are obtained as in Case 2. 



, N9 E 2 El 2c . . . 

y * = ^ + V ~i ~ —sin{2ut + c 2 4.9 



Case 4. The Case of the anisotropic oscillator is well known. The 
trajectories are what is called Lissajous figures. 



2E X 
9^ 



xit) = \l ^sin(3ut + ci) (4.10) 



2E 

y(t) = \l -4sin{ut + c 2 ) (4.11) 



u 2 



Case 5. The potential is 

V = PlyM + (%y/\y\ . (4.12) 
We find the trajectories using the equation (4.4). The bounded trajectories 
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are closed. Examples are shown on Fig 5. 



Case 6. We know from the Section 3 that the only particular case where 
two roots coincide is (3.30). The trajectories for the potentials Vi^x) satisfy 



2E 

y(t) = \l—£*in(ut + c 2 ) (4.13) 



UJ 2 



/dx 
= = t + g (4.14) 
^J-^fx 2 =F ^fxVbTx 2 + 2E 1 

The trajectories were obtained numerically directly from the equations of 
motion. We have closed trajectories and examples are shown on Fig 6. For 
V 3 (x) the potential is a shifted harmonic oscillator, so the trajectories are 
ellipses. 

Case 7. The potential is 

V = a 2 \y\ + b 2 ^/\x~\ (4.15) 



We find the trajectories using equation (4.4). The bounded trajectories are 
closed. They are shown on Fig 7. 

Case 8. There is the particular case 

V = a 2 \y\ + b 2 \x\ (4.16) 

We use eq. (4.4) to calculate the trajectories. Again the bounded 
trajectories are closed. An example is shown on Fig 8. 
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5 Conclusion 



The main results of this article are the study of the algebras of integrals of 
motion, presented in Sections 2 and 3 and the investigation of classical 
trajectories in Section 4. 

In all 8 cases of superintegrable systems, separating in cartesian coordinates 
and allowing a third order integral of motion, the integrals of motion 
generate a cubic Poisson algebra. In many cases this polynomial algebra is 
reducible that is it is a consequence of the existence of a simpler algebraic 
structure. 

In 4 cases the simplest underlying structure is a Lie algebra. Thus in Case 
1, the isotropic harmonic oscillator, the integrals of motion generate the 
u(2) algebra with the Hamiltonian as its central element. The third order 
operator B and the fourth order one C lie in the enveloping algebra of u(2). 
In Case 5 the potential V(x,y) = (3f\/Jx\ + P^yAvl 1S nonanalytical at the 
origin. In this case the Poisson bracket of the second and the third order 
integrals is a constant. Hence A,B and C generate the Heisenberg Lie 
algebra. 

Case 7 with V(x,y) = a 2 \y\ + b 2 -\/\x\ is similar to Case 5 in that A, B and 
C generate the Heisenberg Lie algebra (with C constant). In both of these 
cases the Lie algebra is actually 

{H}Q{A,B,C} (5.1) 

i.e. the direct sum of a one-dimensional Lie algebra, generated by the 
Hamiltonian, and the Heisenberg algebra. 

In Case 8 we have V(x,y)= ay + V(x) with V(x) satisfying the cubic 
equation (3.41), the integrals generate a solvable decomposable Lie algebra 
(3.45). The Hamiltonian H commutes with the elements of the 
two-dimensional solvable algebra. 

The Case 2 and 3 are different. The Hamiltonians (3.6) and (3.12) are 
actually quadratically superintegrable and the third order integral B is the 
Poisson bracket of two second order ones. These second order integrals give 
rise to a quadratic Poisson algebra 15 ' 16 and the cubic algebra of our Section 
3 is an algebraic consequence of the quadratic one. 
Finally, two irreducible cases remain. One is the Case 4, the anisotropic 
harmonic oscillator (3.17). We have a genuine third order polynomial 
algebra (3.19), (3.20) of polynomial integrals H,A,B and C of order 2,2,3 
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and 4, respectively. Somewhat artificially, we can construct the Lie algebra 
u(2) of eq. (3.21), using rational functions of these integrals. 
Case 6 is again irreducible, i.e. we obtain the genuinely cubic polynomial 
algebra (3.33), (3.34), again involving two second order, one third order and 
one fourth order integral. The condition for this algebra to close is precisely 
that V(x) should satisfy the quartic equation (3.28). 

To our knowledge, no classification of polynomial algebras exist, not even of 
quadratic ones, still less of cubic algebras. We can however see that the 
linear parts of the algebras of Case 4 and Case 6 are not isomorphic as Lie 
algebras. That of Case 4 is solvable, that of Case 6 is simple. 
Trajectories are studied in Section 4 with examples on Fig 1,...,8. The most 
important result is that all finite trajectories are closed. We note that the 
potentials of Case 1 and 4 are analytic everywhere, so superintegrability 
implies periodic motion 26 . In cases 2 and 3 the potentials are singular along 
one or both coordinate axes, analytic elsewhere. In Case 5 the potential is 
nonanalytic with nonisolated branchpoints along the axes. That not 
withstanding, bounded trajectories turn out to be closed. 
The functions V\^{x) in Case 6 are nonanalytic for x 2 = —b for b>0. Similar 
statements apply for Case 7. In Case 8 we have studied a psecial case when 
the cubic equation (3.41) has a double root. In this case we either have a 
linear potential V(x,y) = ax + j3y with no bounded trajectories, or the 
nonanalytical one V(x,y) = a 2 \y\ + b 2 \x\ (4.16) that we have investigated. 
Our conclusion is that the bounded trajectories are always closed for the 
superintegrable potentials, wheter they are analytical, or not. 
An investigation of the properties of superintegrable systems with a third 
order integral of motion is in progress. 
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Figure captions 

Figl.A trajectory for V = ^(x 2 + y 2 ). Parameter oj 2 = 2, v xo = —1.5, 
x = 5, vyo = -1.2, y Q = -2, t=[0,400] 

Fig2.A trajectory for ^(x 2 + y 2 ) + ^ + p. Parameter cu 2 = 2,b=2 and 
c=3, v xo = -1.5, x = 5, = -1.2, y a = -2, t=[0,400] 

Fig3.A trajectory for ^-{Ax 2 + y 2 ) + \ + ex. Parameter u 2 = l,b=2 and 
c=3, v xo = -1.5, x = 5, v yo = -1.2, y Q = -2, t=[0,400] 

Fig4.A trajectory for ^(9x 2 + y 2 ). Parameter uj 2 = 1 and d=3, v xo = —1.5, 
x = 5, v yo = -1.2, y = -2, t=[0,400] 

Fig5.A trajectory for V = /3fy/\x\+ Parameter £i = i?2 = = 1 

Fig6.A trajectory for jg (26 + 5x 2 ± 4x^6 + x 2 ). Parameter cu 2 = 1 and 
d=3, u xo = -1.5, x = 5, ^ = -1.2, y Q = -2, t=[0,400] 

Fig7.A trajectory for V = a 2 \y\ + & 2 a/[x[. Parameter Ei = E 2 = k = 1 

Fig8.A trajectory for V = a 2 \y\ + b 2 \x\. Parameter E\ — E 2 — k — 1 
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